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High-resolution,  deep-tow  multichannel  seismic  data  are  used  to  investigate  1ft©  detailed  structure  of  sediments  containing 
methane  hydrate.  These  data  support  thick,  laterally  extensive  layers  of  methane  hydrate-bearing  sediment  unde'  lain  by  a  bottom 
simulating  reflector  (BSR)  and  spatially  discontinuous  zones  of  hydrate  within  the  rndiments  above  the  BSR  den'1’  •~,e  no  BSR 
is  present.  These  data  resolve  normal  faults  which  extend  from  the  surface  throuoh  the  BSR  with  apparent  o?!:>  s  o<  up  to  20  m 
A  phase  inversion  identified  at  the  top  of  the  BSR  shows  that  the  material  immediately  beneath  the  BSR  has  anomalously  low 
velocity,  consistent  with  a  layer  of  sediment  containing  free  methane  gas.  The  fault  offsets  along  the  and  consistent  with  a 
pressure  change  of  -200  kPa  (-2  bars)  across  the  fault  zone. 
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Fine  Structure  of  Methane  Hydrate-Bearing  Sediments  on  the  Blake 
Outer  Ridge  as  Determined  From  Deep-Tow  Multichannel  Seismic  Data 

M.  M.  Rowe  and  J.  F.  Gettrusi 


Naval  Research  Laboratory  Siennis  Space  Center  Mississippi 

High-resulution,  deep-tow  multichannel  seismic  data  are  used  to  investigate  the  deluded  structure  of 
sediments  containing  methane  hydrate  These  data  support  thick,  laterally  extensive  layers  ol  methane 
hydrate-bearing  sediment  underlain  by  a  bottom  simulating  reflector  (BSR)  and  spatially  discontinuous 
zones  of  hydrate  within  the  sediments  above  the  BSR  depth  where  no  BSR  is  present  These  duia 
resolve  normal  faults  which  extend  from  the  surface  through  the  BSR  with  apparent  offsets  of  up  to  20 
m.  A  phase  inversion  identified  at  the  top  of  the  BSR  shows  that  the  material  immediately  beneath  the 
BSR  has  anomalously  low  velocity .  consistent  with  a  layer  of  sediment  containing  free  methane  gas 
The  fault  offsets  along  the  BSR  are  consistent  with  a  pressure  change  of  -  20(1  kHa  t  -2  burs!  across  the 
fault  zone. 


Introduction 

In  1989  the  Naval  Oceanographic  and  Atmospheric  Re¬ 
search  Laboratory,  now  part  of  the  Naval  Research  Labo¬ 
ratory  (NRL).  conducted  a  deep-tow'  multichannel  seismic 
study  on  the  Blake  Outer  Ridge  (Figure  1)  that  included 
profiles  to  resolve  the  detailed  structure  of  methane  hydrates 
which  have  been  identified  tn  that  region  \Shipley  et  al.. 
1979;  Pauli  and  Dillon,  1981;  Markl  and  Bryan.  1983J.  The 
detailed  vertical  and  lateral  distribution  of  gas  hydrates  has 
been  difficult  to  assess  because  the  hydrate  is  stable  only 
within  the  upper  sediments  (to  800  m)  in  relatively  deep 
water  (3000  to  3600  m  in  this  region)  and  are  usually 
identified  indirectly  through  association  with  a  bottom  sim¬ 
ulating  reflector  (BSR)-  The  BSR  is  a  high-amplitude  reflec¬ 
tion  horizon  that  follows  seafloor  topography  and  is  thought 
to  be  the  reflection  from  the  interface  between  hydrate- 
bearing  sediments  at  the  base  of  the  hydrate  stability  zone 
and  gas-saturated  sediments  beneath  them  [Bryan.  1974: 
Shipley  et  al..  1979;  Dillon  et  al..  1980:  Pauli  and  Dillon. 

mi}. 

Gas  hydrate  is  a  specific  type  of  ciathrale.  a  crystalline 
cage  structure,  found  in  deep  ocean  sediments  which  con¬ 
sists  of  a  solid  lattice  of  water  molecules  in  which  each 
“cage"  within  the  lattice  encloses  a  gas  molecule  [Sloan, 
1990J.  Approximately  98%  of  the  gas  making  up  the  hydrate 
in  the  ocean  bottom  is  methane  [Hollister  et  al.,  1972; 
Mclver,  1974;  Kvenvolden  and  McDonald.  1985],  The  meth¬ 
ane  hydrate  buried  in  marine  sediment  contains  high  con¬ 
centrations  of  methane  contained  in  the  hydrate  structure  (as 
much  as  170  m 3  of  methane  gas  in  1  m 3  of  methane  hydrate 
[Baker,  1974;  Kvenvolden  and  Barnard,  1983oJ).  In  addi¬ 
tion,  large  amounts  of  free  (undissolved)  methane  are  be¬ 
lieved  to  exist  trapped  at  the  base  of  the  methane  hydrate 
stability  zone  [Claypool  and  Kaplan.  1974;  Dillon  et  al., 
1980],  This  trapped  methane  may  present  a  significant  haz¬ 
ard  to  drilling  [Miller  et  al..  1991J. 

Methane  hydrate-bearing  sediments  have  distinctive  phys¬ 
ical  properties  which  differentiate  them  from  water-  or 
gas-saturated  sediments.  Compressional  velocities  from  2.0 
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to  3  8  km/s  have  been  measured  in  situ  using  sonobuovs  in 
drill  hole  measurements,  and  in  laboratory  samples  of  meth¬ 
ane  hydrate-saturated  sediment  [Stoll  el  al..  1971;  Brain. 
1974;  Kvenvolden  and  McDonald .  1985:  Mathews  and  con 
Huene.  1985).  Previous  studies  of  methane  hydrate-bearing 
marine  sediments  have  shown  that  the  methane  hydrate  mat 
occur  as  a  solid  mass,  as  a  filling  in  fractures  and  pore 
spaces,  or  as  thin  lenses  within  the  sediment  [Brooks  et  al.. 
1983:  Mathews  and  von  Huene.  1985J.  and  it  may  coexist 
with  methane  dissolved  in  pore  fluids  at  variable  depths 
within  the  methane  hydrate  stability  zone  [Lancelot  and 
£  win  p,  1972J.  The  zone  of  hydrated  sediments  may  be  thin, 
limited  to  the  base  of  the  hydrate  stability  zone  or  may 
extend  from  the  base  of  the  stability  *.one  to  the  lop  of  the 
sediment  column  [MacLeod.  I982J 

The  deep-low  geometry,  source  and  receiver  —35(1  m 
above  the  seafloor  (Figure  2).  and  high-frequency  (250-650 
Hz)  source  signal  used  for  this  study  allow  us  to  resolve 
sediment  structural  details  <5  m  tn  thickness  within  the 
upper  600-700  m  of  sediments.  Lateral  resolution  of  sedi¬ 
ment  structure  of  ~21  m  is  achievable  using  single-channel 
seismic  sections  (i.e..  sections  composed  of  a  single,  com¬ 
mon  offset  trace  from  each  shot  gather).  The  seismic  lines 
recorded  by  NRL's  Deep-Towed  AcoustitVGeophysics 
System  (DTAGS)  are  on  the  western  flank  of  the  Blake 
Outer  Ridge,  at  29°5()'N  latitude.  75°35'W  longitude.  135-165 
kin  west  of  Deep  Sea  Drilling  Project  (DSDPl  sites  102.  103. 
104,  and  533  (Figure  1).  This  region  is  on  the  southwestern 
margin  of  the  area  in  which  a  BSR  has  been  observed  in 
seismic  sections  by  Pauli  and  Dillon  !  1981}  and  Markl  and 
Bryan  [1983], 

The  data  presented  here  show  that  the  compressional 
velocities  within  the  upper  —650  m  of  sediment  are  appro¬ 
priate  for  materials  containing  methane  hydrate,  though 
these  regions  are  not  necessarily  associated  with  a  BSR. 
However,  the  features  diagnostic  of  methane  hydrate¬ 
bearing  sediment  resolved  with  these  seismic  data  occur  in 
laterally  and  vertically  discontinuous  regions  that  appear  to 
be  larger  and  closer  together  where  a  BSR  is  observed 

Deer-Tow  Seismic  Pm* 

DTAGS  multichannel  seismic  system  (Figure  2)  consists 
of  a  Helmholtz  transducer  source  (205  dB  re  1  microPa  (o  l 
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Fig.  I .  Map  showing  relative  locations  of  DTAGS  multichannel  seismic  lines  A  anti  B.  shown  in  Figure  ?.  and  DSDP 
sites.  Water  depth  in  the  urea  within  the  inset  is  3600-3800  m 

m)  and  a  linear  array  containing  24  hydrophone  groups,  array  are  used  to  determine  array  geometry  (see  Gcttrust  rt 

Hydrophone  group  offsets  are  from  137  to  620  m  with  a  «/.  1 19881  and  Gcttrust  and  Rots  1 1  WO]  for  further  discussion 

group  separation  of  21  m.  A  single  calibrated  hydrophone  of  the  system).  Prior  to  the  survey,  conductivity  - 

located  57  m  from  the  source  records  the  source  reference  temperature-depth  (CTD)  measurements  were  made  to  ob- 

s ignat ure.  The  source  signal,  a  125  ms  linear  chirp  from  tain  estimates  ol  the  compressional  velocity  within  the  water 

250-650  Hz.  was  fre  every  30  s,  or  -21m.  The  data  were  column  over  the  full  water  depth,  this  serves  as  "ground 

digitized  at  3125  samples/s  <0.00032  s/sample).  Source  and  truth"  for  water-wave  velocity  estimates  made  with 

array  are  towed  —350  m  above  the  seafloor  in  ocean  depths  DTAGS 

up  to  6000  m.  Depth  sensors  at  the  source  and  along  the  Duta  processing  procedures  are  similar  to  those  applied  to 
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Fig.  2.  Deep-Towed  Acoustics/Geophysics  System  geometry  when  deployed  to  operating  depth.  350-450  m  above  the 
sea  floor.  Digital  data  are  telemetered  to  the  ship  via  the  coaxial  low.  cable 


standard  multichannel  seismic  data;  however,  there  are 
differences  imposed  by  the  deep- tow  techniques  used  to 
collect  these  data.  Because  array  depth  is  not  constrained,  a 
static  correction  is  applied  to  each  channel  to  place  the 
observations  on  a  datum.  As  the  array  tilt  is  always  <  1 1°.  the 
maximum  time  corrections  for  the  array  tilt  are  relatively 
small  (<65  ms). 

The  ship  speed  at  the  slow  ( — 2-0  knots)  tow  speed 
required  by  the  DTAGS  system  is  not  uniform  enough  to 
permit  common  midpoint  (CMP)  processing.  Therefore  all 
data  presented  here  have  been  processed  in  24-fold  common 
shot  point  (CSP)  format.  The  use  of  CSP  gathers,  gathers 
made  up  of  all  24  traces  from  a  single  shot,  will  smooth  the 
velocity-depth  functions  estimated  from  the  data  over  a 
~24Q-m  window  at  the  seafloor.  If  CMP  gathers  had  been 
used,  lateral  resolution  of  the  velocity  structure  would  have 
been  ~21  m  as  each  trace  in  the  gather  would  contain 
reflections  from  the  midpoints  between  the  source  and 
receiver. 

The  data  are  deconvolved  against  the  signal  recorded  at 
the  reference  hydrophone.  It  is  possible  to  use  the  direct 
water  path  signals  because  the  deep-tow  geometry  excludes 
“ghost"  reflections  from  the  water  surface  from  the  data 
window  used  for  these  studies. 

Static  corrections  are  applied  after  signal  deconvolution  to 
remove  the  effects  of  the  array  tilt.  Initial  stacking  velocity 
estimates  are  then  mad*-  using  semblance  analysis.  The 
initial  stacking  velocity,  the  stacking  velocity  of  the  water- 
sediment  interface,  is  the  compressions!  velocity  within  the 
water  column.  This  velocity  estimate  is  compared  with 
(independent)  compressions!  velocity  estimates  based  on 
direct  water-wave  travel  times  to  the  hydrophone  groups  and 
water  column  compressions!  velocity  obtained  from  the 
CTD  measurements.  The  interval  velocity  estimates  lor  the 
water  column  presented  here  are  within  \%  of  those  ob¬ 
tained  from  the  concurrent  direct-path  measurements  and 
CTD  measurements. 

Stacking  velocity  estimates  are  made  using  at  least  every 
third  shot;  therefore  with  a  (nominal)  21 -m  shot  separation, 
estimates  are  made  at  63-m  intervals  along  the  profile.  Even 


with  this  relatively  short  spatial  separation,  there  are  places 
where  lateral  heterogeneity  requires  even  finer  sampling.  In 
such  regions,  velocity  estimates  are  made  using  near-trace 
(offsets  from  137  to  368  ml  gathers  from  consecutive  shot 
gathers  (21  m  apart),  decreasing  the  amount  of  lateral 
smoothing  to  - 100  m.  Velocity  estimates  are  checked  using 
depth  sections,  created  by  transforming  two-way  travel  time 
to  depth  using  these  velocity  estimates.  This  is  a  valuable 
test,  as  abrupt  changes  in  horizon  depths,  inconsistent  with 
the  sediment  structure  shown  in  the  single-channel  time 
section,  will  occur  where  velocity  estimates  are  poor.  Depth 
sections  also  allow  us  to  compare  our  estimated  depth  to  a 
major  reflection  horizon  with  any  coring  te.g..  DSDP)  data 
taken  near  our  profiles. 

Observations 

Two  multichannel  seismic  lines  are  presented  in  this 
paper.  The  first,  line  A.  the  eastern  line  in  Figure  1 ,  is  2.5  km 
long  and  trends  from  north-south.  The  second,  line  B.  the 
western  line  in  Figure  I ,  is  7  .9  km  long  and  trends  from  west 
to  east.  These  lines  are  located  approximately  150  km  west 
of  DSDP  sites  102.  104.  and  533.  on  the  southwestern  flank 
of  the  Blake  Outer  Ridge.  The  seismic  data  are  characterized 
by  numerous  closely  spaced  horizontal  reflectors  extending 
more  than  10  km  across  the  two  sections  (Figure  3).  Similar 
closely  spaced  reflectors  are  observed  in  conventional  mul¬ 
tichannel  seismic  data  collected  in  this  area  ! Hollister  cl  at.. 
1972.  Villon  pi  nl ..  1976;  Shipley  cl  a!..  1979;  Mark I  and 
Bryan.  1983).  The  deep-tow  data  re»_al  that  these  closely 
spaced  reflectors  are  broken  by  high-angle  normal  faults  less 
than  I  km  apart  extending  to  the  maximum  depth  of  the 
section,  approximately  0.65  s  below  the  water-sediment 
interlace  (Figures  3  and  4i. 

Cores  recovered  at  DSDP  sites  102.  103.  104.  and  *'33  on 
the  Blake  Outer  Ridge  (Figure  1 1  consist  of  uniform  hemipe- 
lagic  silty  clay  interspersed  with  thin  beds  of  silty  to  sandy 
material  to  at  least  660  m  depth  [Hollister  c r  a!..  1972; 
Sheridan  et  at..  19831.  This  stratigraphy  is  the  resuli  of 
fluctuations  in  the  location  and  strength  of  the  contour 
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Fig.  3.  («)  Neur-truce  single  channel  seismic  section  from  line  A.  characterized  by  thm.  closely  spaced  reflectors 

and  downs  lope  normal  faulting.  A  high-amplitude.  shallow  reflecting  horizon  is  noted  by  the  arrow  on  the  left-hand  side 
The  arrow  on  the  right-hand  side  shows  the  BSR  (hi  Ne.ir-trace  single  channel  seismic  section  from  line  B  The 
laminated  reflecting  horizons  in  this  section  correlated  well  with  the  reflectors  in  line  A  though  the  lines  are  5  km  apart 
and  orthogonal  to  each  other.  (Frequency  of  the  reflection  data  has  been  reduced  for  display .  > 


currents  which  controlled  deposition  of  these  late  Miocene 
to  Holocene  sediments  on  the  flanks  of  the  Blake  Outer 
Ridge  [£ii’ing  and  Hollister.  1972:  Mark!  and  Bryan.  19831. 
Common  offset  DTAGS  data  resolve  the  small-scale  (-5  m 
vertical.  -21  m  horizontal)  changes  in  the  sediment  struc¬ 
ture  and  reflectivity  [.BoiWe.v  et  al..  1991  ].  The  closely  spaced 
acoustic  horizons  seen  in  Figure  3  are  reflections  from  the 
small  erosional  unconformities  [Mark!  and  Bryan.  19831  and 
minor  changes  in  sediment  composition  and  density  which 
are  apparent  in  the  core  samples. 

The  water  depth  where  these  deep-tow  seismic  data  were 
recorded  ranges  from  3660  to  3680  m:  this  is  within  the  range 
of  water  depths  at  DSDP  sites  102,  103.  and  104  to  the  east. 
As  water  depth,  geologic  environment,  and  sediment  type 
are  similar  at  these  DSDP  sites  and  within  the  DTAGS 
survey  area  (Figure  I).  the  maximum  depth  of  methane 
hydrate  stability  at  the  deep-tow  seismic  site  will  be  within  a 
few  tens  of  meters  of  the  maximum  depth  of  hydrate  stability 
at  DSDP  sites  102  and  104:  that  is.  615-620  m  | Hollister  et 
al..  1972). 

The  strong,  continuous  BSR  at  1.28  to  1.30  s  (all  travel 
times  presented  here  are  two-way  travel  times).  0.64  to  0.65 
s  below  the  walcr-sediment  interface  (Figure  3 a)  appears 
only  in  seismic  line  A  (Figure  I ).  This  reflecting  horizon  has 
the  high-amplitude  characteristic  of  a  BSR  and  occurs  at  the 
same  arrival  time  as  the  BSR  observed  in  seismic  sections 
recorded  near  DSDP  sites  102  and  104.  0.6  s  J Lancelot  and 
Ewin/i.  t972|.  We  show'  below-  that  the  depth  of  this  BSR  is 
—650  m  below  the  water-sediment  interface. 

In  seismic  line  B  there  is  no  discernable  BSR-type  reflect¬ 
ing  horizon  prescr.:  :n  the  DTAGS  multichannel  seismic 
data,  though  seismic  horizons  above  the  BSR  observed  in 
the  DTAGS  data  car.  be  correlated  exactly  troin  line  A.  the 
eastern  section,  through  line  B.  the  western  section  (Figure 
3).  A  low-amplitude  reflection  of  limited  extent  may  he 
discerned  between  3.5  and  4.5  km.  0.65  s  below  the  water- 
sediment  interface  in  Figu.e  3 b.  coincident  with  the  arrival 
time  of  the  BSR  in  Figure  3</. 

Comparison  of  Figures  3u  and  3 h  shows  that  a  shallow, 
high-amplitude  acoustic  reflector  ( indicated  by  the  arrow  on 


the  left  side  of  Figure  3 a  a:  0.83  si  on  line  A  marks  the  top 
of  a  zone  of  decreased  reflectivity  relative  to  the  overlying 
reflectors.  This  zone  of  decreased  reflectivity  extends  down 
to  the  BSR.  Reflectivity  does  not  appear  to  decrease  below 
this  horizon  on  line  B 

On  line  A,  a  series  of  parallel  normal  faults  dipping 
between  44'  and  60“  extend  from  the  top  of  the  sediment 
column  down  through  the  BSR.  the  amount  of  upparent 
offset  increasing  w-ith  increasing  depth  (Figure  4)  Fault 
offsets  in  the  BSR  have  not  been  resolved  by  surface-tow 
seismic  data.  The  high  resolution  achievable  with  DTAGS 
enables  the  offsets  in  the  BSR  to  be  clearly  imaged  in  the 
seismic  data. 

Progressively  larger  offset  with  depth  indicates  fault  slip 
along  these  zones  of  weakness  has  been  occurring  contem¬ 
poraneously  with  sediment  deposition  since  late  Miocene  or 
early  Pliocene  time.  The  sense  of  the  faulting,  downslope 
along  the  ridge  flank,  is  consistent  with  gravity-driven  sliding 
of  the  ridge  flank  sediments.  Growth  faults  like  these  have 
been  observed  on  the  continental  slope  to  the  north  [ Car¬ 
penter.  1981:  Dillon  et  al..  I983J.  Slope  movement  along 
some  of  the  growth  faults  in  these  regions  is  attributed  to 
failure  of  the  sediments  along  a  layer  of  weak,  gas-rich 
sediments  underlying  the  more  rigid  methane  hydrate  layer 
1  Carpenter.  1981). 

These  faulls  are  observed  in  line  E  (Figure  3b).  though 
they  do  not  exhibit  similar,  parallel  offset  since  this  seismic 
line  is  close  to  parallel  to  the  strike  of  the  fault  planes.  The 
strike  of  the  fault  planes  of  these  gravity-driven  growth  faults 
is  nearly  parallel  to  the  seafloor  contours. 

S>T"  'u  Vixotiw  Asm  ys*;; 

Compressional  velocity  estimates  obtained  from  the  east¬ 
ern  and  western  DTAGS  seismic  lines  show  vertical  vari¬ 
ability  as  small  as  a  few  tens  of  meters  in  extent  and  lateral 
variability  of  the  order  of  a  hundred  meters  (Plate  I ).  In  both 
sections,  compressional  velocity  varies  little  within  the 
upper  0.1 8  s  below  the  water-sediment  interface,  increasing 
only  slightly  from  the  compressional  velocity  at  the  sediment 
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Fig.  4.  Near-trace  single  channel  seismic  section  shown  in  Figure  3 a  expanded  to  shim  parallel  normal  faults  ulipping 
44°  to  60°)  penetrating  the  BSR.  (Frequency  of  the  reflection  data  has  been  reduced  for  display .  I 


surface,  from  1.53  to  1.65  km/s.  Impedance  contrasts  in 
these  sediments  are  driven  by  density  changes  within  thin 
(~  I  wavelength)  reflecting  layers  associated  with  smali-scale 
erosional  unconformities  [ Markl  and  Bryan,  1983).  Large, 
abrupt  changes  in  compressional  velocity  across  reflecting 
horizons  are  not  consistent  with  either  our  compressional 
velocity  estimates  or  DSDP  velocity  data. 

At  0.80  to  0.83  s  on  line  A  (the  top  of  the  zone  of  low 
reflectivity)  the  compressional  velocity  increases  sharply 
(Plate  I  a);  velocity  increases  from  —1.65  km/s  to  >2.00 
km/s  over  0.09  s.  Between  this  marker  horizon  (0.20  s  below 
the  seafloor)  and  the  BSR  (0.60  s  below  the  sediment 
surface),  the  compressional  velocity  increases  to  2.4  to  2.7 
km/s  (Plate  la).  Constant  velocity  normal  move-out  (NMO) 
corrected  shot  gathers  (Figure  5)  demonstrate  the  increasing 
compressional  velocity  with  depth.  On  !me  B  this  steep 
increase  in  velocity  is  not  observed  continuously  ucross  the 
section  though  locally,  in  regions  of  the  order  of  a  hundred 
meters  in  extent,  there  is  a  steep  compressional  velocity 
gradient  in  which  the  velocity  increases  to  2.3  km/s  (Plate 
1  b).  These  high-velocity  regions  occur  within  the  same  depth 
range  as  the  zones  of  low  reflectivity  in  line  A  (Figure  3«). 


The  BSR  al  the  base  of  line  A  (Figure  3«)  is  a  wide 
(—0.005  to  —0.010  s)  reflecting  horizon,  indicating  that  it  is 
composed  of  reflections  from  the  top  and  bottom  of  a  layer, 
having  a  thickness  of  the  order  of  one  to  two  wavelengths. 
Estimates  of  the  compressional  velocity  w>ithin  this  thin 
layer  cannot  be  made  since  a  change  in  NMO  corrections 
between  1.0  and  3.0  km/s  at  these  depths  and  velocities 
would  be  of  the  order  of  a  single  data  sample  interval  A 
distinct  phase  inversion  is  observed  at  the  top  of  the  BSR.  at 
1.292  s  but  not  at  the  lower  interface  al  1.300  s  (  Figure  h«i, 
indicating  that  the  compressional  velocity  of  the  material 
immediately  beneath  the  upper  interface  of  the  BSR  is 
significantly  lower  than  the  compressional  velocity  of  the 
overlying  sediment. 

Figure  6 h  shows  a  synthetic  seismogram  consisting  of  a 
shot  gather  containing  traces  at  the  same  offsets  as  show  n  in 
Figure  6 a.  137-242  m.  The  best  fit  model  was  a  low-velocity 
(1.3  km/s)  layer  4-m-thick  overlying  higher- velocity  (2.1 
km's)  material  (Figure  6r>.  Similar  estimates  of  the  velocity 
and  thickness  of  this  thin  layer  were  obtained  by  Miller  a  al. 
1 1991 1  for  the  BSR  beneath  the  inner  slope  of  the  Pern-Chile 
trench . 
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Plate  1.  to)  Two-dimensional  interval  velocity  model  determined  from  the  multichannel  data  shown  in  Figure  3u. 
The  velocity  model  extends  down  to  the  BSR.  at  ~  1.3  s.  ( b )  Two-dimensional  interval  velocity  model  derived  from  the 
multichannel  seismic  data  shown  in  Figure  ib.  Below  1.24  s  there  are  no  strong,  continuous  reflecting  horizons  for 
making  velocity  estimates. 


Lateral  variability  in  compressional  velocity  in  line  B  is  as 
high  as  0.3  km/s  over  distances  of  the  order  of  100  tn,  the 
minimum  lateral  smoothing  window  achievable  with  these 
data,  (Plate  1  b).  Isovelocity  regions  as  narrow  as  -150  m 
wide  exist  which  the  sound  velocity  remains  relatively  low. 
1. 5-1. 7  km/s.  from  the  sediment  surface  to  the  bottom  of  the 
section  Examples  of  such  zones  are  between  0.9  and  1.05 
km  and  between  2.8  and  3.2  km  in  Plate  1  h.  These  zones  are 
bounded  by  faults  in  the  seismic  sections. 

Figures  la  and  lb  are  depth  sections  derived  from  the 
data  shown  in  Figures  3 a  and  3 b,  converting  from  time  to 
depth  using  the  interval  velocity  functions  shown  in  Plates 
1  a  and  1  b .  The  closely  spaced  reflecting  norizons  above  the 
BSR  can  be  traced  from  line  A  through  line  B.  even  though 
the  lines  are  orthogonal  and  there  is  considerable  lateral 
variability  in  the  velocity  functions.  Diffractions  from  the 


sharp  edges  of  the  faulted  horizons  have  not  been  migrated. 
The  BSR  occurs  approximately  650  m  below  the  seafloor  and 
varies  in  thickness  between  approximately  3.5  and  7.0  m 
(assuming  the  compressional  velocity  within  this  low  veloc¬ 
ity  layer  is  1.3  km/s).  This  is  consistent  with  the  BSR  depth 
of  615-620  m  estimated  from  drilling  results  from  DSDP  sites 
102  and  104  | Hollister  el  a!..  1972). 

Spatially  averaged  compressional  velocity-depth  profiles 
derived  from  our  data  (Figure  8«)  resolve  vertical  variability 
in  compressional  velocity  on  a  scale  of  a  few  tens  of  meters. 
These  fine-scale  velocity  changes  are  not  resolved  b\  Hamil¬ 
ton’s  (1979)  empirical  compressional  velocity-depth  function 
(Figure  8 a)  for  the  upper  600  m  of  fine-grained,  terrigenous, 
seafloor  sediments.  Hamilton's  [  1979)  velocity  function  pre¬ 
dicts  a  much  higher  compressional  velocity  gradient  (>T$) 
in  the  upper  100  m  ot  sediment  than  is  obtained  from  DTAGS 
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Fig.  5.  Sample  shot  gather  from  line  A  {0.95  km.  Figure  So),  normul  move-out  (NMO)  corrected  using  the 
isovelocity  functions  given  above  each  gather.  A  given  horizon  will  be  horizontal  across  the  gather  when  NMO 
corrected  using  the  correct  stacking  velocity.  Horizons  below  1 .08  s  flatten  at  significantly  higher  velocities  (1.71  to  1 .8 1 
km/s)  than  horizons  above  1 .00  s.  This  shows  that  a  steep  velocity  gradient  exists  between  1 .00  and  1 .08  s.  (Frequency 
of  the  reflection  data  has  been  reduced  for  display.) 

data  (close  to  0/s),  even  in  sediments  which  do  not  contain 
methane  hydrate  [Gettrust  et  ul.,  1988],  The  vertical  vari¬ 
ability  in  compressional  velocity  as  a  function  of  depth 
shown  by  our  data  is  not  resolved  by  the  velocity  profiles 
derived  from  conventional,  surface-towed,  multichannel 
seismic  data  such  as  the  data  from  the  Blake  Outer  Ridge 
discussed  by  Dillon  anti  Pauli  [1983]  (Figure  8 h).  This 
variability  is  not  shown  by  the  compressional  velocity  profile 
estimated  by  relating  drilling  breaks  at  DSDP  sites  102  and 
104  to  reflecting  horizons  in  surface-towed  single-channel 
data  [Hollister  et  al. ,  1972]  which  is  smoothed  over  hundreds 
of  meters  in  depth  (Figure  8 h). 

Methane  Hydrate  Occurrence 

The  two-way  travel  time  to  the  BSR  in  the  DTAGS  data. 

O.fi  s  below  the  sediment  surface,  coincides  with  the  two- 
way  travel  time  to  the  BSR  observed  in  surface-towed 
seismic  data  by  Hollister  et  ul.  [  1972),  Tucholke  et  al.  1 1977). 
and  Shipley  et  al.  [  1979).  The  subbottom  depth  of  the  BSR  in 
Figure  8c  correlates  with  the  maximum  depth  of  methane 
hydrate  stability  predicted  from  the  phase  relation  of  Kven- 
volden  and  Barnard  [  1983c]  and  MacLeod  [1982],  This 
depth  also  corresponds  to  the  depths  at  which  drilling  breaks 
occurred  at  DSDP  sites  102  and  104. 

The  geothermal  gradient  estimated  from  the  depth  of  the 


BSR  in  the  DTAGS  data,  using  the  relation  of  Tucholke  et  al. 
[1977]  and  sediment  density  of  1.5  g'em '  ( Paulas .  1972]  is 
approximately  3.7°C/100  m.  This  estimate  is  close  to  the 
3.6°C/100  m  measured  value  reported  by  Kvenvolden  and 
Barnard  ]  1 9836]  at  DSDP  site  533  and  is  w'ithin  the  3,5(T  to 
3.95°C/!00  m  range  given  by  Stall  and  Brvan  11979]  as 
reasonable  estimates  of  the  geothermal  gradiem  in  sediment 
containing  methane  hydrate. 

The  consistently  high  (>2.0  km/s)  compressional  velocity, 
the  decreased  reflectivity  (used  as  a  diagnostic  for  estimating 
hydrate  layer  thickness  by  Dillon  et  al.  1 1991  ]l  between  250 
and  050  m  depth,  and  the  continuous  BSR  indicate  that 
methane  hydrate  forms  a  thick,  continuous  layer  within  the 
sediment  between  approximately  250  and  650  nt  depth  below 
the  seafloor  in  the  area  covered  by  line  A. 

Velocity  analysis  of  the  data  from  the  western  Blake  Outer 
Ridge  line,  line  B.  shows  evidence  of  localized,  discontinu¬ 
ous  methane  hydrate  formation  even  though  a  distinct  BSR 
is  not  present  in  the  seismic  sections.  The  homogeneous 
lithology  indicated  by  the  core  samples  from  DSDP  sites  102 
and  104  j Hollister  et  al..  1972]  is  inconsistent  with  the 
laterally  highly  variable  compressional  velocity  and  reflec¬ 
tivity.  Lateral  variability  in  velocity  observed  in  sediments 
in  which  methane  hydrate  is  no!  present,  such  as  the 
Bermuda  Rise,  is  significant!!  lower  than  that  observed  here 
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Fig.  6.  (a)  Near-trace  shot  gather  at  1.0  km  in  line  A  showing  the  phase  inversion  at  the  top  of  the  BSR  1 1  2V  s) 

lb)  Near-trace  synthetic  seismogram  which  hesl  fit  the  gather  in  Figure  6«.  This  seismogram  indicates  that  ’he  phase 
is  not  inverted  below  the  BSR.  at  1.30  s.  (c)  Best  fit  compressional  velocity  model  used  to  generate  the  synthetic 
seismogram  in  Figure  fib. 


[ Gettrust  et  aL,  1988).  The  high  lateral  variability  in  sedi¬ 
ment  compressional  velocity  observed  within  the  predicted 
zone  of  hydrate  stability  (extending  to  650  m  subbottom 
depth)  in  line  B.  1.7  to  2.1  km/s  over  distances  of  the  order 
of  100  m,  is  evidence  that  regions  of  methane  hydrate¬ 
bearing  sediments  as  limited  as  100  m  in  extent  coexist  with 


similarly  sized  regions  of  sediments  containing  pore  fluid  and 
dissolved  methane  (if  methane  is  preseni)  bui  no  hydrate. 
The  existence  of  such  regions  has  been  proposed  b\  Lance¬ 
lot  and  Ewing  |1972)  and  Galimov  and  Kvcnvoiden  [1983). 

The  localized,  discontinuous  development  of  methane 
hydrate  within  the  sediment  may  be  due  to  lateral  variation 


Fig.  7.  (a)  Depth  section  from  line  A  (Figure  3«)  converted  to  depth  using  the  velocity  function  shown  in  Plate  l«. 

(b)  Depth  section  from  line  B  (Figure  3fe)  converted  to  depth  using  the  velocity  function  shown  in  Piute  lb.  ( Frequency 
of  the  reflection  data  has  been  reduced  for  display.) 
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Fig.  8.  (a)  Spatially  averaged  compressions!  velocity-depth  functions  from  lines  A  and  B  as  compared  to 

Hamilton's  1 1979]  empirical  compressional  velocity-depth  curve.  Velocity  function  for  line  A  helms  the  BSR  depth  tf>sti 
m)  is  based  on  modeling  results,  lb)  Compressional  velocity-depth  functions  obtained  from  surface-lowed  multichannel 
seismic  data  \Dillnn  and  Pauli,  1983]  and  DSDP  site  102  and  104  drilling  rate  data  referenced  to  reflecting  horizons  in 
single-channel  seismic  data  \Hollixter  et  al..  1972). 


in  methane  concentration  as  the  gas  is  generated  from 
organic  material  in  situ  {Brooks  et  a!..  1983:  Dillon  and 
Pauli.  1983  j  or  to  loss  of  methane  from  regions  of  the 
sediment  which  do  not  presently  contain  hydrate.  Destabili¬ 
zation  of  the  methane  hydrate  along  the  fault  zones  may 
occur  due  to  the  drop  in  pressure  and  increase  in  frictional 
temperature  during  rupture.  As  a  result,  large  quantities  of 
methane  may  be  released  from  the  sediments  near  the  fault 
zone  over  a  short  period  of  time  [  Reed  et  al..  1 990| .  Between 
rupture  events  these  faults  are  barriers  to  lateral  (and  verti¬ 
cal)  gas  migration,  preventing  hydrate  from  reforming  in  the 
methane-depleted  regions  within  the  zone  of  hydrate  stabil¬ 
ity. 

Absence  of  a  high-amplitude,  continuous  BSR  in  line  B  is 
interpreted  to  be  the  result  of  an  overall  decrease  in  the 
sediment  methane  concentration  between  the  lines  A  and  B. 
The  sporadic  occurrence  of  methane  hydrate  inferred  from 
the  velocity  function  shown  in  Plate  I  b  and  the  seismic 
sections  from  line  B  (Figures  3 b  and  lb)  may  be  character¬ 
istic  of  the  transition  from  a  region  in  which  methane  hydrate 
is  present  in  a  thick,  continuous  layer  marked  by  a  BSR  at 
the  base,  as  seen  in  Figure  3«.  to  an  area  where  less  methane 
is  present  and  discontinuous  regions  containing  high  concen¬ 
trations  of  methane  form  localized  pockets  of  methane 
hydrate-bearing  sediment. 

N  ature  or  the  BSR 

Interpretation  of  the  BSR  as  the  equilibrium  phase  bound¬ 
ary  between  methane  hydrate  and  methane  gas  and  water 
appears  to  conflict  with  the  displacement  of  the  BSR  shown 
in  the  seismic  section  in  Figure  3.  Geochemical  evidence 
suggests  that  the  BSR  could  be  caused  by  reflection  from  a 
dense  ankerite  or  siderite  horizon,  precipitated  as  a  result  of 
the  high  concentrations  of  methane  in  the  pore  fluid  beneath 
the  methane  hydrate  phase  boundary'  [ Lancelot  and  Ewing. 
1972;  Ctaypnol  and  Kaplan.  1974:  Matxumoto.  1989:  Esikov 
and  Paxhkina.  1990).  Thus  while  related  to  the  location  of 


the  maximum  depth  of  stability  of  methane  hydrate,  the 
actual  reflecting  horizon  is  caused  by  the  impedance  contrast 
between  the  dense  carbonate  formed  in  situ,  such  as  the 
"very  hard  layer  or  nodule"  of  ankerite  dribed  al  the  depth 
of  the  BSR  at  DSDP  site  104  {Hollister  et  al..  I972|.  and  the 
overlying  gassy  hemipelagic  marine  sediment.  This  carbon¬ 
ate  layer  would  be  offset  by  faulting  but  would  transgress 
time  horizons  within  the  sediment  in  the  same  way  as  the 
methane-hydrate  phase  boundary 

However,  the  high-amplitude  phase  inversion  at  the  BSR 
(Figure  ba)  is  evidence  that  the  impedance  difference  is 
caused  by  a  significant  decrease  in  compressional  velocity  at 
the  BSR  rather  than  a  velocity  increase.  This  behavior  is 
consistent  with  a  model  in  which  the  BSR  is  the  reflection 
from  the  top  of  a  lever  of  sediment  containing  methane  gas 
rather  than  the  reflection  from  a  hard  layer  of  carbonate.  The 
coda  of  the  BSR  reflection  is  best  modeled  (Figure  bh)  by  a 
thin  (—4  m)  layer  of  low  compressional  velocity  (1.3  km's) 
sediment  at  the  depth  of  the  BSR  overlying  sediment  with  a 
significantly  higher  velocity  (2. 1  km/s)  (Figure  be).  Thicken¬ 
ing  of  the  BSR  updip  along  the  faulted  sections  (Figure  4) 
may  be  due  to  migration  of  the  free  gas  upward  along  the 
base  of  the  impermeable  hydrate. 

The  thin,  gas-rich  layer  below  the  BSR.  suggested  by  the 
model,  would  have  low  shear  strength,  providing  a  surface 
along  which  slope  failure  may  occur.  Growth  faults,  like 
those  seen  in  line  A  (Figure  4).  are  characteristic  of  slope 
failure  along  a  deeper,  less  rigid  sediment  layer.  Dillon  et  al. 
f  1 983 ]  observed  such  faulting  in  w'hich  the  less  rigid  layer 
was  salt.  Carpenter  [1981]  observed  these  grow'th  faults  in 
which  failure  occurred  along  a  layer  of  gas-rich  sediment. 
Dissociation  of  the  hydrate  al  the  base  of  the  BSR.  due  to 
small  pressure  or  temperature  changes,  may  facilitate  slip 
along  this  zone  of  weakness  [  Pauli  et  al.,  1991 ;  Dillon,  1991  ]. 

Distinct  offsets  in  the  equilibrium  phase  boundary  marked 
by  the  BSR  are  consistent  w’ith  the  discontinuity  in  the 
compressional  stress  field  across  the  fault  zone  |  Hafner, 
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1951),  The  increase  in  depth  to  the  BSR  suggests  pressure  is 
higher  at  the  BSR  depth  on  the  downthrown  side  of  the  fault 
relative  to  the  upthrown  side.  Higher  pressure  on  the 
downthrown  sides  of  these  faults  is  consistent  with  the 
observed  displacement  along  these  growth  faults  {Hapicr. 
1951).  The  pressure  increase  across  the  fault  from  the 
upthrown  side  to  the  downthrown  side  max  abruptly  in¬ 
crease  the  maximum  depth  of  methane  hydrate  stability  as 
much  as  ~  20  m  at  the  fault  zone.  This  pressure  effect  would 
decrease  with  increasing  distance  from  the  fault,  causing  the 
depth  to  the  BSR  in  this  case  to  decrease  to  the  south 
between  the  faults  (Figure  4). 

Given  the  water  depth  ( — - 37<Kf  m)  and  depth  to  the  BSR 
(—650  m)  in  the  area  studied,  the  phase  diagram  of  kven- 
\ widen  and  Barnard  |  1 98 3 n .  Figure  1 1.  indicates  that  a 
change  in  pressure  across  a  fault  of  -  200  kPa  (-2  bars) 
would  be  consistent  with  a  BSR  displacement  of  15-20  m. 
This  differential  stress  is  within  the  bounds  of  the  shear 
strength  of  saturated  clavs  at  these  depths  [(ran;.  1985.  p 
135). 

CONCI  l.'SIONS 

Evidence  from  deep-towed  multichannel  seismic  data 
(compressional  velocity  structure,  variability  in  reflectivity, 
phase  inversions  at  the  RSR)  shows  that  methane  hy'dratc 
may  occur  continuously  within  a  -400-m-thick  layer  of 
sediment  extending  down  to  the  maximum  depth  of  hydrate 
stability,  characterized  by  low  reflectivity  and  underlain  by  a 
BSR.  or  it  may  he  limited  to  regions  on  the  order  of  a  100  m 
in  extent  where  high  concentrations  of  methane  exist  in 
discontinuous  patches  Patchy  development  of  methane  hy¬ 
drate  may  he  due  to  loss  of  methane  from  the  sediment  or  to 
a  decrease  in  the  amount  of  methane  generated  in  situ. 

Velocities  of  the  order  of  2.5  kmN  may  be  accounted  for 
by  lithologic  changes  within  the  sediment  such  as  the  devel¬ 
opment  of  dense  carbonate  or  chert  layers  However,  the 
uniform  lithology  of  the  Blake  Outer  Ridge  sediments  be¬ 
tween  the  water-sediment  interface  and  depths  exceeding 
ftflO  m.  as  shown  by  the  cores  from  DSDP  sites  102.  103.  and 
104.  and  the  continuous  nature  of  the  reflecting  horizons 
extending  8  km  across  line  A  to  line  B.  5  km  away,  do  not 
support  lithologic  changes  as  the  cause  of  the  anomalous 
high  velocities  in  either  the  eastern  or  the  western  sections. 

These  high-resolution  data  show  that  the  BSR  has  finite 
thickness  of  5-10  m.  The  phase  inversion  at  the  top  of  the 
BSR  shows  that  this  5-  to  10-m-thick  layer  is  a  low-velocity 
zone  (most  likely  free  gas  hearing)  which  also  is  a  zone  of 
weakness  that  may  he  causally  related  to  the  norma)  faults 
observed  penetrating  the  BSR.  7’he  ~  20-m  fault  offset  re¬ 
solved  with  these  data  along  the  BSR  is  not  resolved  in 
surface-towed  seismic  data.  As  noted  previously,  this  offset 
is  consistent  with  a  reasonable  estimate  of  the  pressure 
differential  across  the  fault.  If  there  is  surficiai  expression  of 
these  faults  that  can  be  resolved  with  side  scan  sonar  images, 
side  scan  sonar  mapping  of  methane  hydrate  regions  could 
serve  as  a  rapid  means  of  refining  estimates  of  the  areal 
extent  of  these  zones:  this,  in  combination  with  additional 
high-resolution  multichannel  seismic  data,  would  provide 
additional  insight  into  the  volume  of  methane  hydrate  exist¬ 
ing  in  the  marine  environment .  Although  only  two- 
dimensional  profiles  were  taken  during  the  expenmcnl  dis¬ 
cussed  here,  the  ability  to  resolve  methane  hydrate  zones 


within  deep-ocean  sediments  ias  demonstrated  with  these 
data)  indicates  thai  it  will  be  possible  to  use  high-resolution 
multichannel  seismic  data  in  the  fuiurc  to  obtain  viable  in 
situ  estimates  of  the  volume  of  methane  hy  drate  and  hence 
the  importance  of  these  features  m  terms  of  energv  sources 
and  their  potential  impact  on  the  greenhouse  effect  ( Paid!  rt 
at..  1991) 
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